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Рассмотрена осесимметричная задача о распространении фронта фазового превращения в твердой порис- 
той среде в условиях граничного теплообмена при наличии теплозащитного слоя. Проведен параметриче- 
СКИЙ анализ, позволяющий определять влияние величины теплового потока на границе, а также теплофи- 
зических свойств изоляционного слоя на показатели волны фазового превращения. 
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Введение. Рассматриваемая задача актуальна в связи со строительством и эксплуатацией зда- 
ний и сооружений в районах Крайнего Севера в условиях вечной мерзлоты, когда обусловлен- 
ные природными и искусственными факторами процессы плавления и кристаллизации содер- 
жащейся в грунте влаги приводят к нарушению устойчивости объектов вплоть до аварийных 
ситуаций. В частности, это имеет место при бурении и эксплуатации нефтяных и газовых сква- 
жин, когда прилегающий к устью скважины и содержащий включения влаги горный массив под- 
вержен влиянию потока тепла со стороны самой скважины. Наличие кольцеобразного теплоза- 
щитного слоя, примыкающего к внутренней границе среды, оказывает влияние на параметры 
волны фазового превращения влаги в ней. Причем здесь имеет значение не только толщина 
этой прослойки, но и теплофизические свойства наполняющего её вещества. Схематично задача 
представлена на рис.1. 
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Рис.1. Схема зон: 1 — внутренняя цилиндрическая стенка; 2 — кольцевой теплозащитный 
слой; 3 — граница между теплозащитным слоем и пористой средой; 4 — зона пребывания 


влаги в твердом виде; 5 — зона пребывания влаги в жидком виде; И, - радиус внутрен- 


ней цилиндрической полости: Й -— толщина теплозащитного кольцевого слоя: И ф те 


кущий размер зоны фазовых превращений 
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Математическая модель задачи основана на тех же, что и в [1], представлениях. Твер- 
дая пористая среда содержит равномерно распределенные включения влаги, которые в рас- 
сматриваемом температурном диапазоне могут находиться в одном из двух агрегатных состоя- 
ний: либо в твердом, либо в жидком. Соответствующие зоны твердой пористой среды, обла- 
дающие при осреднении по объему различными теплофизическими свойствами, будем иденти- 
фицировать как зоны твердых и жидких включений. Считается, что на макроскопическом уров- 
не фазовое превращение происходит на поверхности раздела этих зон. Теплофизические свой- 
ства вещества каждой из зон, включая зону теплового барьера 2 (см. рис.1), полагаются одно- 
родными по пространству и неизменными во времени. 

Математическая модель задачи в безразмерных переменных и величинах включает в се- 
бя уравнения распространения тепла для каждой из трех зон 
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той 7 в момент времени 1; х, — радиус внутренней цилиндрической полости; с, ^,, р -— 


удельные теплоемкости, коэффициенты теплопроводности и плотности веществ соответ- 
ственно, индекс «с» относится к теплоизоляционному слою, индекс «м» относится к зоне 
пребывания влаги в твердом состоянии, «т» - к зоне, где влага находится в жидком со- 


СТОЯНИИ, й — толщина теплозащитного слоя, отделяющего содержащий В СОВОКУПНОСТИ 


влагу массив от кольцевой полости. 
Указанным выше областям соответствуют следующие интервалы пространственной пе- 


ременной 7: х <г<к-+й - теплозащитный слой; и, +й <Т < - слой, подверженный фазо- 

вым превращениям (далее Р-слой); к+й<и<и + й + ’, - зона жидких включений; 

№) + А +7, <г < - зона твердых включений. Здесь А =7, +й + ’, — координата фронта фазово- 
ох, 

го превращения (радиус зоны кристаллизации влаги); А=——, а — коэффициент теплоотдачи 


[@ 


от хладагента, циркулирующего вдоль цилиндрической стенки к твердому массиву; 
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Т., Т - масштабные значения температуры, первое из которых соответствует характерной 


температуре потока хладагента, второе — температуре твердой среды на значительном удале- 
нии от центральной оси. Следует заметить, что безразмерный комплекс 4 связан прямо пропор- 
циональной зависимостью с числом Нуссельта указанного выше омывающего стенку потока. 

Задача дополняется условиями сопряжения на границах, разделяющих зоны. На непод- 
вижной границе «теплозащитный слой — Р-слой» должно выполняться условие непрерывности 
потока тепла: 
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где х.„„ — момент времени, начиная с которого в твердой пористой среде будет распростра- 
ФП 


няться волна фазового превращения. На свободной границе раздела фаз выполняются 
условия непрерывности температуры и теплового баланса: 
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где У =—*^; © =-^ —, а=——" ,Т (ТГ <Т, <Т ) - температура фазового перехода; 
д тт’ 1 = —т в (1. <, <Т.) ратура ф р 


4, — теплота фазового перехода в расчете на единицу массы вещества среды; Г - ско- 





рость перемещения фронта фазового превращения. Очевидно, что х,(т)=0 на проме- 


жутке времени те [0;т.„ |. 
Как видно, математическая модель (1)-(7) включает в себя параметры теплозащитного 
слоя й, Л., С. , которые наряду с параметрами А-слоя С,, Л, определяют динамику фронта 


м 
фазового превращения. 

Численный анализ задачи. Для численного анализа, как и в [1], с помощью разбиения фи- 
зического пространства на контрольные объёмы, выделения узловых точек и интегрирования 
уравнений (1) по каждому из контрольных объёмов строится дискретный сеточный аналог. Сле- 
дует заметить, что разбиение производится так, чтобы одна из узловых точек располагалась в 
точности на границе между теплозащитным слоем и Е-слоем. 

В качестве базового набора параметров были приняты следующие значения входящих в 
задачу констант: /,, =1,25; Л, =1; с, =0,74; с,=1; 4,=5,08; й=0,2; Л, =0,76; с, =0,74; 
Т, =0,75; А=10. 

На приводимых ниже рисунках представлены результаты численного анализа рассмат- 
риваемой задачи, позволяющего оценить влияние параметров теплоизоляционного слоя на 
процесс распространения зоны фазовых превращений в твердой пористой среде. Так, на рис.2 
представлены результаты расчетов, показывающие влияние коэффициента теплопроводности 
^., вещества, составляющего теплозащитный слой, на радиальную глубину проникновения зоны 


т 


фазового равновесия в область, занятую веществом с жидкими включениями. Как и следовало 
ожидать, с уменьшением А не только снижается темп нарастания слоя с влагой в твердом со- 


стоянии, но и увеличивается время задержки х„„ процесса распространения волны фазовых 
превращении. 
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Рис.2. Динамика зоны фазового превращения при различных значениях коэффициента 
теплопроводности Л. теплозащитного слоя 


На рис.3 представлена зависимость безразмерного радиуса зоны кристаллизации влаги 
в среде от величины коэффициента теплопроводности теплозащитного слоя за определенный 
промежуток времени, который в рассматриваемом случае имеет значение 1,3. Из рисунка вид- 
но, что существует значение /^., =)., , ниже которого в течение указанного промежутка времени 


фазовые превращения в пористой среде не наблюдаются. 
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Теплопроводность /, 
Рис.3. Зависимость размера слоя с кристаллизовавшейся влагой от значения коэффициента 
теплопроводности вещества теплозащитного слоя за время {=1,3 
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Разработанный алгоритм численного решения позволяет производить количественную 
оценку влияния теплопроводящих свойств изоляционного материала при фиксированной тол- 
щине защитного слоя и, в том числе, определять пороговое значение ).., ниже которого рас- 


пространение процесса фазового превращения за указанный промежуток времени не происхо- 
дит. 

Влияние режима теплообмена на цилиндрической стенке осевой полости на динамику 
процесса фазовых превращений показано на рис.4. Режим конвективного теплообмена на гра- 
нице между потоком хладагента (например, газовым потоком) и кольцевым теплозащитным 
слоем определяется величиной коэффициента теплоотдачи А, который тем больше, чем выше 
интенсивность циркуляции хладагента вдоль цилиндрической стенки. Анализ результатов чис- 
ленных расчетов показывает, что при фиксированных прочих параметрах задачи, начиная с не- 
которого значения, дальнейшее повышение коэффициента теплоотдачи (соответственно, уве- 
личение скорости обтекания стенки) не оказывает существенного влияния на темп распростра- 
нения процесса фазового перехода в пласте. Это обусловлено тем, что для больших А время 
охлаждения границы «хладагент - твердая среда» до температуры омывающего полость потока 
невелико и динамика процесса практически определяется скоростью теплопереноса в твердой 
фазе. 
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Рис.4. Динамика зоны фазового превращения при различных значениях коэффициента 
теплоотдачи А на цилиндрической стенке 


Еще более отчетливо подтверждает данное обстоятельство рис.5, на котором приводит- 
ся зависимость радиуса зоны отвердевания влаги от величины теплоотвода за фиксированное 
время (1=1,3). В конкретном случае эти зависимости позволят обосновать режим циркуляции 
хладагента вдоль границы твердой пористой среды. 
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Рис.5. Положение фронта фазового равновесия по истечении времени {=1,3 
в зависимости от значения коэффициента теплоотдачи на граничной цилиндрической стенке 
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Толщина теплозащитного слоя й 
Рис.6. Расчетное положение фронта фазового равновесия и толщина зоны с твердыми включениями 
на момент времени {=1,3 в зависимости от толщины теплозащитного слоя 


На рис.б представлены результаты расчетов, позволяющие оценить влияние толщины 
теплозащитного слоя на глубину зоны кристаллизации. Как видно, с увеличением толщины изо- 
ляционной рубашки поперечный размер образовавшегося за фиксированное время (1=1,3) слоя 
с твердыми включениями (кривые 1? и 2?) монотонно падает. В то же время зависимость от той 
же величины радиуса фронта кристаллизации (кривые 1° и 2°) содержащейся в среде влаги, 
включающего совокупную толщину защитного и «мерзлого» слоев, имеет немонотонный харак- 
тер. При этом наблюдающийся в этой зависимости минимум сдвигается в сторону больших зна- 
чений толщины защитного слоя с уменьшением коэффициента теплопроводности /,, последне- 


го. С практической точки зрения, выявленные и имеющие количественное описание закономер- 
ности могут служить основанием для подбора теплоизоляционной прослойки, обладающей оп- 
тимальными по отношению к выдвигаемым критериям параметрами. 
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Рис.7. Время распространения фронта кристаллизации на расстояние Аи = 0,2 


в зависимости от величины теплообмена на цилиндрической стенке: 1 -— И=0,2; 2 — Й=0,1 


Время промерзания (1$ 





Теплопроводность /, 
Рис.8. Время распространения фронта кристаллизации на расстояние Л» = 0,2 в зависимости 


от значения коэффициента теплопроводности Л. теплоизоляционного слоя: 1 — й=0,2; 2 — й=0,1 


На рис.7 и 8 представлены результаты расчета времени {, распространения процесса 
радиальной кристаллизации влаги в среде на определенное расстояние Ди (в рассматриваемом 
случае Л, = 0,2.) для двух значений толщины й тепловой изоляции. Рис.7 иллюстрирует зави- 


симость этого времени от интенсивности теплообмена на внутренней цилиндрической стенке, а 
рис. 8 — влияние теплопроводности защитного слоя на продолжительность процесса распро- 
странения фронта кристаллизации на указанную длину. 
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Заключение. Анализ полученных результатов показывает, что для управления процессом рас- 
пространения зоны фазового равновесия претерпевающего фазовый переход компонента твер- 
дой пористой среды существует несколько возможностей. В частности, добиться требуемого те- 
плового состояния твердого массива можно с помощью выбора подходящего режима конвек- 
тивного теплообмена на его границе, используя наряду с этим изоляционный материал с опре- 
деленными теплофизическими и структурными свойствами. 
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ЕКОМТ 1М НЕТЕКОСЕМЕОЦ$ ОО РОКОЦ$ МЕОТОМ 1М ТЕКМ$ ОЕ ВОЧМОАКУ 
НЕАТ ТКАМ$ЗЕЕВ АМО РКЕЗЕМСЕ ОЕ СЕМТКАЕ НЕАТРКООЕ КТМС ГАУЕК 


А.М. СОГОМММА, Т.$. КАЕУАКТМ 
(УКЖа Зае Тесптса! Упмег$Ку), 
В.А. МЕУООВЕ 

(Роп Чае Тесптса! Упмег$Ку), 


\.М. РОЗНКТМ 


(УКЖа Зае Тесптса! Упмег$Ку) 


Ап ах5уттес ргоМет оп ргораданоп о! пе рвВазе КапзоптаНоп топ т те 5о/а рогоиз$ тей ит 
п сопаоп$ ор пе ВБоипаагу ВеаЁ гап$Гег п {Ве ргезепсе оЁ! Пеаргоо! [ауег 15 сопэ!аегей. Тре рагате- 
[С апау$/5 (паЁ реттй5 {о деегтте {пе тйиепсе о! {Пе ПеаЁ Пои/ гаёе оп {Пе Боипоаагу, апа те {пегторрВуз!- 
са/ ргорегиех о! {Пе тпзи/аЕптд Гауег оп {Ве рагатеег$ о! {пе рпПа5е гапГогтайоп игауе [5 тасе. 

Кеуигога$: 50! рогоиз теит, рразе гапГогтаНоп, опЁ ргорадаНоп, ахуттеётс етап ргоМет, ВеаЕ- 
ргоо! /ауег, питегса! о/иНоп. 
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